Finite temperature behaviour of glueballs in Lattice Gauge Theories 
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We propose a new method to compute glueball masses in finite temperature Lattice Gauge The- 
ories which at low temperature is fully compatible with the known zero temperature results and as 
the temperature increases leads to a glueball spectrum which vanishes at the deconfinement transi- 
tion. We show that this definition is consistent with the Isgur-Paton model and with the expected 
contribution of the glueball spectrum to various thermodynamic quantities at finite temperature. 
We test our proposal with a set of high precision numerical simulations in the 3d gauge Ising model 
and find a good agreement with our predictions. 
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While the physics of glueballs in pure LGTs at zero 
temperature is by now rather well understood [TH3] a 
similar level of understanding for the finite tempera- 
ture behaviour of the glueball spectrum is still lacking. 
The standard method used to compute finite T glueball 
masses [HE] is to measure the correlator of the glueball 
operator (for instance a simple plaquette, if one is inter- 
ested in the ++ glueball). along the compactified time 
direction of length i . The glueball spectrum obtained in 
this way turns out to be almost constant as the tempera- 
ture increases and seems not to be affected by the decon- 
finement transition: glueball masses were measured even 
deeply in the dcconfmed phase showing values similar to 
the zero temperature ones (or slightly smaller, depend- 
ing on the procedure adopted in the calculation) [3J [5] . 
However this picture is unsatisfactory for at least two 
reasons. 

First, one of the most successful phenomenological 
descriptions of glueballs is the well known Isgur-Paton 
model [6\ This model and its recent generalizations [7 
is able not only to predict the general structure of the 
spectrum (like, for instance, the fact that the mass of the 
2 ++ state is lower than the mass of the 1 ++ ) but also its 
fine details and shows a remarkable agreement with the 
lattice estimates. In this model glueballs are considered 
as " rings of glue" , kept together by the same string ten- 
sion which appears in the interquark potential and should 
thus vanish at the deconfinement point when the string 
tension vanishes. If we trust this picture, then it could 
be used also to predict the T dependence of the glueball 
spectrum for low temperatures. In fact the Isgur Paton 
model predicts values of the zero temperature glueball 
masses m,i(T = 0) as adimensional ratios Tii(0)/-\/cr(0) 
(where <r(0) is the zero temperature string tension). In a 
finite temperature setting we expect the same ratios, but 
with er(0) substituted by the finite temperature string 
tension u{T) i.e.: 
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dependence proposed in [H [5] for this range of tempera- 
tures. 

Second, recently, very precise estimates of various ther- 
modynamic quantities have been obtained both below 
and above T c . in pure lattice gauge theories in d = 3 + 1 
[Sl[9] and in d = 2 + 1 QUI [IT] dimensions. For T < T c the 
thermodynamics of these theories is very well described 
in terms of a gas of glueballs which are the only degrees 
of freedom of the theory in this regime. For T > T c 
the thermodynamics is well described by a gas of free 
gluons (and accordingly the thermodynamic observables 
scale as N 2 ). If glueballs were present also in the decon- 
fined phase they would give an additional contribution to 
the thermodynamic observables, leading to results fully 
incompatible with the lattice measurements. 

These observations suggest that with the current 
method to extract finite T glueball masses [H [5] one is 
probably measuring some other finite size scale of the 
model whose relation with the glueball spectrum is sim- 
ilar to the relation which exists between the spacelike 
string tension a s and the finite temperature string ten- 
sion cr(T). Indeed also <r s , which is extracted from space- 
like Wilson loops, is almost constant for T <T C , increases 
for T > T c and it is well known to be completely unre- 
lated to the finite temperature string tension a(T) which 
is instead extracted from Polyakov loop correlators. 

In this letter we propose an alternative prescription to 
evaluate finite T glueball masses, compatible with the 
above observations. The most direct way to ensure the 
expected finite T behaviour is to construct an observable 
with the correct quantum numbers so as to be coupled 
in the continuum limit to the glueball states, using only 
Polyakov loops so as to have the correct dependence on 
the finite temperature string tension. The simplest pro- 
posal is to choose a pair PP^ of nearby Polyakov loops 



M(x) = P(x)P i (x + a) 



(2) 



which is a decreasing function of T. This expectation 
is in complete disagreement with the almost constant T 



where a denotes the lattice spacing and P(x) is the 
Polyakov loop located at the space point x. Then the 
glueball mass will be extracted looking at the large R 
behaviour of the connected correlator of two M(x) oper- 
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FIG. 2. Effective string description of our proposal. 



FIG. 1. The glueball correlator discussed in the text. 



ators as depicted in figjl] 

G(R, T) = (M(O)M(R)) - (M) 2 fl "„, e - m "^ R (3) 

The space-like version (usually denoted as " torelon pair" ) 
was used in [21 [3] as part of the operator basis to obtain 
the T — glueball spectrum. In the proposed interpre- 
tation, this set up is new. 

Let us also stress that this is not the only possible 
choice, for instance in non-abelian gauge theories an 
equivalent interesting possibility would be the Wilson 
loop obtained joining the two Polyakov loops with two 
space like links at t = and t = N t (this choice is obvi- 
ously equivalent to our proposal for abelian LGTs). 

The nice feature of our proposal is that it has a nat- 
ural interpretation in terms of the effective string model 
of pure gauge theories. It is the four point correlator of 
four closed effective strings (see fig. [2]). The external 
legs correspond to the four Polyakov loops (which are 
described as closed strings due to the compactification 
of the time direction), while the glueballs are the excita- 
tions of the closed string joining together the four legs. 
As mentioned in the introduction this proposal is strongly 
based on our intuition of the glueball dynamics coming 
from the Isgur-Paton model. It might be useful to make 
more explicit this connection. If we could perform a sec- 
tion in the four strings function as depicted in fig. 2, the 
effective string description of the flux distribution within 
the section would be given by the Isgur Paton model. 
Accordingly we expect that all the glueballs (indepen- 
dently of their quantum numbers) would flow within the 
closed string, of which they would represent different ra- 
dial or rotational excitations. In the large R limit only 
the lowest mass survives, but in principle, looking at the 
subleading exponentials for lower values of R one could 
recover also the remaining states of the spectrum. While 
the radius of the four external legs is fixed to be the in- 
verse of the finite temperature, the radius of the internal 
closed string coincides with the glueball radius tq which 
is one of the parameters of the Isgur-Paton model. 

Despite all these interesting features there is appar- 
ently a major problem with this proposal. In fact, due 
to dimensional reduction, one expects that any mass 
scale extracted from an observable of this type should 
scale in the vicinity of the deconfinement transition as 



m s (T) ~ o(T)/T which, as it is easy to see, is pretty 
different from the expected scaling behaviour of eq.Q. 

We shall see below in a concrete example how this 
problem can be addressed. Indeed, as we shall see, for 
any T < T c , a (glueball) mass with the correct scaling be- 
haviour is always present in the spectrum of G(R, T), but 
as the deconfinement transition is approached this mass 
becomes subleading and the large R behaviour is dom- 
inated by a different mass scale (whose physical mean- 
ing we shall discuss below) with the "wrong" scaling be- 
haviour m s (T) ~ a(T)/T. However, even in the vicinity 
of T c , there is always a suitable range of values of R in 
which the glueball mass, even if subleading, can be un- 
ambiguously observed. 



Test in the 3d gauge Ising model 

In order to test our proposal we computed the mass of 
the lightest glueball in the 3d gauge Ising model. This 
choice has two relevant advantages. 

First, very precise estimates exist for the zero tem- 
perature spectrum 12J with which we can compare our 
results in the low T limit. In particular we know that 
in the range of f3 values that we study in this paper 
too =3.15(5)^0(0) [Hj. 

Second, using dimensional reduction |13j and the fact 
that the 3d gauge Ising model has a second order de- 
confinement transition in the same universality class of 
the 2d Ising magnetization transition, we can predict the 
behaviour of the correlator G(R, T) in the vicinity of T c 
using results borrowed from the exact solution of the 2d 
Ising model. In this limit our observable becomes equiv- 
alent to the energy-energy correlator in the high temper- 
ature phase of the 2d Ising spin model. From the exact 
solution of the model we know that the large R behaviour 
of the function should be dominated by a new mass scale 
m s which is exactly twice the fundamental mass of the 
model, which, from dimensional reduction, is known to 
be <j(T)/T (see for instance the discussion in sect. 2. 3 of 
|14)h Thus, as anticipated, we expect in this limit 



m s {T) ~ 2a(T)/T 



(4) 



We performed three sets of simulations at different val- 
ues of the gauge coupling (corresponding to ~- = 5.67a, 
8a and 12a respectively [H]) in order to test scaling cor- 
rections. For each value of (3 we chose a value of the lat- 
tice size in the spatial direction L s large enough to make 



finite size effects negligible, and studied various values 
of compactified time direction N t = 1/T in the range 
T < T c . For each value of T we evaluated the correlator 
G(R, T) for several values of R. We also evaluated for 
each T in a separate simulation the finite temperature 
string tension a(T) (using the methods discussed in [T5] ) 
so as to be able to construct the scaling functions eqs. (fTl) 
and Q . A few details on the simulations are reported in 
tab. m 
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i 


L s 


N t 


R 


0.743543 


5.67 a 


90 


7,8,9 


6<R<20 


0.751805 


8 a 


90 


9,10,11,12,13,14,20,56,64 


8<R<22 


0.756427 


12 a 


120 


20 


12 < R < 33 



TABLE I. For each of the three j3 values we report the cor- 
responding critical temperature T c and the values of L s , Nt 
and R that we studied. 
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FIG. 3. 



and 
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plotted as a function of S- for /3i = 



la{T) ^/<x(T) 

0.743543, p 2 = 0.751805 and /? 3 = 0.756427. The two curves 
correspond to the two expected scaling behaviours: mo(T) ~ 
3.15y/(r(T) and m s (T) = 2a(T)/T. 



We found two different behaviours. For low values of 
T (in our simulations the threshold was T < 0.6T C ) the 
data were perfectly fitted by the following expression 



G(R,T) = a (T)- 



-m (T)R 



R 



(5) 



with good x 2 values in the whole range of values of R 
that we considered. The data were so precise that we 
could also confirm the presence of the expected 1/y/R 
prefactor. The values of TOo(T) extracted from the fits 
are reported in the last few lines (for each j3) of tabJTT] 
and in fig [3] and turned out to follow exactly the expected 
behaviour eq.(nj), with a value of the glueball mass in 
good agreement with the T = value m ~ 3.15i/cx(0) 
obtained in [T2"] . 

For high values of T (i.e. in our case T > 0.6T C ) it 
turned out to be impossible to fit the data using eq.((5]). 
Reasonable x 2 values could only be obtained discarding 
the low R values of the correlators and using a different 
fitting function: 



G(R,T)=a s (T)- 



,(T)R 



R 2 



(6) 



We use the notation m s to stress the fact that this mass 
was obtained using a different fitting function. The pref- 
actor 1/R 2 is exactly what one would expect for the 
energy-energy correlator in the 2d Ising model and the 
mass m s extracted form this fit scales exactly as sug- 
gested by eq.Q. In full agreement with the expectation 
of dimensional reduction, not only the T dependence but 
also the fact that m s is exactly twice the fundamental 
mass of the model is perfectly reproduced by the data 
(see tab. [IT] and fig'J3]). 

This explains the large R behaviour of G(R,T), how- 
ever in order to include also the small R data in the fit 
it turned out to be mandatory to use a two exponentials 



fitting function: 



G(R,T) = a (T)- 



-m {T)R 



R 



o.(T)- 



-m s (T)R, 

~R^ 



(7) 



It is well known that this type of fits is very delicate. 
In our case we used the following procedure: we kept the 
lowest mass m s fixed to the value obtained in the large R 
limit and fitted again the data keeping as free parameters 
only ao (T) , a s (T) and mo (T) we could in this way fit all 
the data with a reduced x 2 of order unity. It is important 
to stress that the identification of the subleading expo- 
nential was facilitated by the very different behaviour of 
the two prefactors, by the wide range of values of R that 
we used in the fit and by the fact that these data were 
not cross-correlated since, due to the algorithm that we 
used (see [15]). each value of R was obtained in an inde- 
pendent simulation. These observations should be taken 
into account when trying to reproduce our results in other 
LGTs. The results of our fits are reported in tab. [II] and 
fig[3] The subleading mass Too in the T > 0.6T C region 
turned out to be the natural continuation of the glueball 
mass that we had identified for T < 0.6T C . As it is easy 
to see looking at fig[3] mo follows eq. ([lj up to the highest 
temperatures that we studied with a value m ~ 3.15 in 
agreement with the T = estimate. 

It would be important to understand the effective 
string interpretation of this new scale m s . We have no 
rigorous proof but it is likely that the crossover that 
we observe between too and m s is due to a competition 
between the two minimal surfaces bounded by the four 
Polyakov loops which are compatible with the topology 
of the lattice and of our observable. As the temperature 
increases it becomes less and less costly for the flux tube 
to wind around the periodic boundary conditions lead- 
ing to a minimal surface composed by two parallel flux 
tubes as depicted in fig[4] . This crossover is controlled 
by the competition of two scales: the compactification 




FIG. 4. Minimal surface associated to the m s mass. 



radius 1/T and the glueball radius r . As T increases 
also ro increases (since it is due to the flux tube width 
which is known to increases with T) thus it will certainly 
exist a crossover value of T above which ro > 1/T which 
in the 3d gauge Ising model that we studied turns out to 
be around T c /2. 






a(T) 



,{T)T 



r(T) 



TTTT 



n^JTT 



0.743543 
0.743543 
0.743543 
0.751805 
0.751805 
0.751805 
0.751805 
0.751805 
0.751805 
0.751805 
0.751805 
0.751805 
0.756427 



0.8 

0.7 
0.62 
0.89 

0.8 
0.73 
0.67 
0.62 
0.57 

0.4 
0.14 
0.125 

0.6 



0.00961 
0.01315 
0.01542 
0.00268 
0.00444 
0.00566 
0.00654 
0.00720 
0.00771 
0.00922 
0.01037 
0.01040 
0.00326 



2.03(4) 
1.89(5) 
1.98(5) 
2.2(2) 
1.97(5) 
1.86(8) 



1.39(3) 
1.73(5) 
2.21(6) 
1.02(8) 
1.31(3) 
1.54(7) 



3.1(2) 
3.1(1) 
3.14(8) 
3.0(3) 
2.9(3) 
2.9(1) 
3.3(1) 
3.23(3) 
3.29(5) 
3.25(4) 
3.14(3) 
3.21(2) 
3.29(6) 



TABLE II. Values of a(T), 



m s (T) 
a[T) 



^£1 and ^±. 



Concluding Remarks 

The main message of our analysis is that the mass of 
the lowest glueball (and thus likely the whole glueball 
spectrum) scales at finite temperature as y/a(T) and 
thus is a decreasing function of T/T c and vanishes at 
T = T c . Our results also suggest that the Isgur-Paton 
model is valid also at finite temperature and that its pre- 
dicted scaling behaviour (mo ~ y/a(T)) can be concili- 
ated with the different scaling behaviour predicted by the 
Svetitsky-Yaffe conjecture [T3] thanks to the appearance 
of a new mass scale (m s ~ a(T)/T) which in the vicinity 
of T c dominates the large R behaviour of the correlator. 
This is likely to be a general mechanism. For instance, a 
similar phenomenon was also observed a few years ago in 
the finite T behaviour of the monopole spectrum of the 
random percolation gauge theory (see fig.2 of ref.[l7]l. 
It is also interesting to notice that this new mass scale 
strongly resembles the "spurious states" observed in [5J|3] 
which, in fact, were characterized by a large overlap with 
the torelon pair states (we thank B. Lucini for this ob- 



servation). 

It would be interesting to understand the physical mean- 
ing of to s . Preliminary simulations show that at high 
enough temperatures the picture we have discussed holds 
almost unchanged even if we increase the distance be- 
tween the two nearby Polyakov loops up to a few lattice 
spacings. In this limit our observable describes the inter- 
action of two mesons and, according to the effective string 
picture discussed above, the mass scale m s should mea- 
sure the attractive interaction between quarks and anti- 
quarks belonging to different mesons. Our results show 
that this interaction becomes the dominant contribution 
in the meson-meson correlator as T c is approached from 
below. This agrees with the intuitive picture of decon- 
finement as a "melting" of mesons into individual quarks. 
In our framework this melting transition would be driven 
by the interaction mediated by the mass scale m s . 
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